





ACT-05/97 

CTP-TAMU-10/97 

D OE/ER/40717-40 

hep-ph/970235q 



New supersymmetric contributions to t ^ cV 



^ ; Jorge L. Lopez^, D.V. Nanopoulos^'^, and Raghavan Rangarajan^ 

^ [ ^ Bonner Nuclear Lab, Department of Physics, Rice University 

^ : 6100 Main Street, Houston, TX 77005, USA 

'^ ■ ^Center for Theoretical Physics, Department of Physics, Texas A&M University 

^ '■ College Station, TX 77843-4242, USA 

O . ^Astroparticle Physics Group, Houston Advanced Research Center (HARC) 
The Mitchell Campus, The Woodlands, TX 77381, USA 



^ . Abstract 

We calculate the electroweak-like one-loop supersymmetric contributions to 
k> , the rare and flavor- violating decay of the top quark into a charm quark and a 

;h I gauge boson: t — > cV, with V = j,Z, g. We consider loops of both charginos 

and down-like squarks (where we identify and correct an error in the literature) 
and neutralinos and up-like squarks (which have not been calculated before). 
We also account for left-right and generational squark mixing. Our numerical 
results indicate that supersymmetric contributions to t ^ cV can be upto 5 
orders of magnitude larger than their Standard Model counterparts. However, 
they still fall short of the sensitivity expected at the next-generation top-quark 
factories. 
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1 Introduction 

The discovery of the top quark by the CDF and DO Collaborations []I]| at Fermilab, 
and its subsequent mass determination {rrit = 175 ± 6 GeV) have initiated a new 
era in particle spectroscopy. However, unlike the lighter quarks, the top quark is not 
expected to form any bound states, and therefore its mass and decay branching ratios 
may be determined more precisely, both theoretically and experimentally. Upcoming 
(Run II - Main Injector) and proposed (Run III - TeV33) runs at the Tevatron will 
yield large numbers of top quarks, as will be the case at the LHC, turning these 
machines effectively into 'top factories.' Even though higher precision in the deter- 
mination of rrit is expected (it is already known to 3%), more valuable information 
should come from the precise determination of its branching fractions into tree-level 
and rare (and perhaps even 'forbidden') decay channels. 

The purpose of this paper is to study one class of such rare decay modes: 
t -^ cV, with V = 'y,Z,g. The particular case oi t ^ eg has received some phe- 
nomenological attention recently as a means to probe the scale at which such new 
and unspecified interactions might turn on [Q. Our purpose here is to consider an 
explicit realization of this coupling within the framework of low-energy supersymme- 
try. This is different in spirit from the line of work in Ref. 0, as the effective mass 
scale at which such vertices 'turn on' is determined here by the interactions of pre- 
sumably rather light sparticles. Within supersymmetry, the t -^ cV vertex was first 
contemplated in Ref. 0], where the one- loop QCD-like (loops of gluinos and squarks) 
and electroweak-like (loops of charginos and down-like squarks) contributions were 
calculated. The QCD-like supersymmetric corrections were subsequently re-evaluated 
and generalized in Ref. [Q, which pointed out an inconsistency in the corresponding 
results of Ref. [^. Here we study the electroweak-like supersymmetric contributions 
to t ^ cV. We reconsider the chargino-down-like-squark loops and point out and 
correct an inconsistency, essentially a lack of gauge invariance because of the apparent 
omission of a term, in the corresponding results of Ref. 0. We also consider for the 
first time the neutralino-up-like-squark loops, and include the effects of left-right and 
generational squark mixing. 

Our numerical results indicate that for typical values of the parameters one 
gets a large enhancement over Standard Model predictions of top-quark decays to 
gauge bosons . For the most optimistic values of the parameters the enhancement 
can be as large as five orders of magnitude. However, even for the most optimistic 
values of the parameters, such rare decay channels fall short of the expected sensitivity 
of the next-generation top factories. 

2 Analytical Results 

In this section we obtain the one-loop electroweak-like supersymmetric effective top- 
quark-charm-quark-gauge-boson vertex by considering loops involving charginos and 
neutralinos, including the effects from left-right and generational squark mixing. We 



then present the decay rates of the top-quark to the charm-quark and a gauge boson. 
The invariant amphtude for top-quark decay to a charm-quark and a gauge 
boson can be written as 

M = Mo + 5M (1) 

where Mq is the tree-level amplitude and 5M is the first-order supersymmetric cor- 
rection. As there are no explicit flavor-violating tcV couplings in the Lagrangian 
Mq = 0, whereas 6M is given by 

tSM = u{p2)V^u{pi)e,{k,X), (2) 

where pi,p2, and k are the momenta of the incoming top-quark, outgoing charm- 
quark, and outgoing gauge boson respectively, and e^(fc. A) is the polarization vector 
for the outgoing gauge boson. The vertices V^ may be written as 

y^(tcZ) = -iYiPLFzi + PRF'z^) + k,a>^''{PRFz2 + PLF'z2), (3) 

V^itcj) = -i^^{PLF^i + PRF'^^) + k,a^'''{PRF^2 + PLFl^2), (4) 

V^{tcg) = -tT'^^^{PLFg^ + PRF'g,)+T'^k,a^''{PRF,2 + PLF'g2). (5) 

where as usual we have defined Prl = i(l ± 75) and a^'^ = |[7'^,7'^]. T" are the 
generators of SU{3)c- The form factors -Fi_2 and F^^ encode the loop functions and 
depend on the various masses in the theory. The Feynman rules used to obtain them 
are given in Refs. P, |^ and the corresponding Feynman diagrams are shown in Fig. [l|. 
The vertices are derived assuming pi — p2 — k = 0. 

The form factors for the electroweak-like corrections due to loops involving 
charginos and strange- and bottom-squarks are given by 
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In these expressions i = Z,'y,g, the sums over j,k = 1,2 run over the two chargino 
mass eigenstates, p,e = 2,3 represent strange- and bottom-squarks (we ignore the 
mixing with the down-squark) , and l,7n = 1,2 represent a sum over squark mass 
eigenstates which are obtained from the ql^r gauge eigenstates via: qpi = cos 9 pQpL + 
svo-OpfipR and qp2 = — sin 9pqpL + cosOpCipR. The various B and c functions in the 
above expressions are the well documented Passarino-Veltman functions ||^ (adapted 
to our metric where pj = 7nf)\ the arguments of the B and c functions are indicated 
by the superscripts on the braces in Eqs. (§|)-(0). For example, the arguments of the 
c functions that appear within the first brace of Eq. (|^) are (— pi, fc, rn +,7nq^^, rriq^^) 
while the arguments of the B functions that appear in the third brace of Eq. (^ are 
(— Pi,m ±,mg J. (Note that the Passarino-Veltman functions depend only on the 
square of their arguments.) The Passarino-Veltman functions contain infinities which 
cancel each other out, as they should since there is no t-c-V vertex in the Lagrangian. 
The coefficient functions are given by 
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The chargino mixing matrices Uij and Vij and the generational mixing matrices K, r2 
and B2 which appear in these expressions are defined in Ref. 0, M^ is the diagonal 
matrix {mdjmsyTnf,) and Tql is the squark mixing matrix defined in Ref. |^. 

In deriving the above form factors, we have used the following relations: 
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(The first two equalities in Eq. (^) are only true when one takes the absolute square 
of both sides of the equation and uses Eq. (|29|) .) 

Our results above for the chargino-squark loops disagree with those of Ref. 
in the limit of mr = 0. We have an additional term 
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in F^^. This term is required by gauge invariance, i.e., it is needed to ensure that the 
coefficient of 7^ in Eqs. (H) and and (|^) vanishes for the massless gauge bosons. 

The form factors for the electroweak-like corrections due to loops involving 
neutralinos and top- and charm-squarks are given by Eqs. (|^)-(P) with m^± replaced 
by m^o and with the coefficients A^-Fc replaced by 
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The neutralino mixing matrices Nij and N^j and the generational mixing matrices 
Fi and Bi, are defined in Ref. |^, and M„ is the diagonal matrix {mu,mc,mt). In 
Eqs. (^-(^ the sums over j, k now run from 1-4 over the four neutralino mass 
eigenstates, and p,e,p = 2,3 represent charm- and top-squarks (we ignore the mixing 
with the up-squark); /,m = 1,2 represent a sum over squark mass eigenstates, as 
earlier. The coefficients G and H are unaltered. 

While we have included the charm-quark mass in the form factors above for 
completeness, we set ttIc = hereafter. The supersymmetric electroweak-like contri- 
bution to the decay rates is then given in terms of the form factors obtained above, 
as follows 
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In these expressions each form factor receives contributions from both char gino- squark 
and neutralino-squark form factors. It is not hard to verify that for ttIc = 0, the 
chargino contributions to F[ 2 vanish. (In Ref. [Q there is a factor of vr missing in the 
expression for T{t ^ cZ) and a factor of | missing in r(t — ;> eg).) 



3 Numerical Results 

Before we attempt to evaluate the rather lengthy expressions given above, we would 
like to consider qualitatively the possibility of dynamical enhancements of the loop 
amplitudes. Experience with similar diagrams contributing to the self-energy of the 
top quark in supersymmetric theories ||10| indicates possible large corrections when 



the mass of the top quark equals the sum of the masses of the other particles leaving 
the vertex involving the top quark.Q In the present case we have vertices with top 
quarks and: (i) gluinos and top-squarks in the case of QCD-like contributions (cal- 
culated in Refs. [0, Q); (ii) charginos and down-hke squarks in the case of 'charged' 
electroweak-like corrections; and (iii) neutralinos and up-like squarks in the case of 
'neutral' electroweak-like corrections. Given the presently-known lower bounds on the 
squark (excluding t) and gluino masses {i.e., rriq, rrig > 175 GeV; rriq ^ rrig > 230 GeV 
|]TI|), this type of enhancement might only be present in the third type of contribution 
when nit ~ "^x + ^™ti ' "which requires a light top-squark whose mass is constrained ex- 



perimentally to ?7ij^ > 60 GeV |T^. The 'neutral' electroweak-like contributions might 
also be enhanced by large GIM-violating top-squark-charm-squark mass splittings. 
This latter enhancement is also present in the 'charged' electroweak-like contributions. 
However, the 'charged' contributions fall short of the 'neutral' electroweak-like correc- 
tions for the gluon and photon cases. Interestingly, the 'charged' contribution for the 
Z is higher than in the 'neutral' case. We first address the neutral electroweak-like 
contributions and comment on the charged contributions afterwards. 

The neutral electroweak-like contributions might be enhanced as discussed 
above, but this is subject to other mixing factors in the An — Fn coupling functions 
in Eqs. (p^-(^) being unsuppressed. At the root of this question is whether the 
quark-squark-neutralino couplings might be flavor non-diagonal as a result of their 
evolution from the unification scale down to the electroweak scale. This question 
might be explored by considering the squared squark mass matrices at the electroweak 
scale that are obtained by renormalization group evolution of a universal scalar mass 
at the unification scale j^, |T^: 

X^ = M^/xSj + /xSx,X| (z = 1,2) (50) 

Xl^ = M^ fi^^ll + fi^^lx.Xl + fiflx^Xl (51) 

^2L = M^f^fll + ix^^lx^Xl + fxflx^Xl (52) 

where i = 1 (i = 2) corresponds to up-type (down-type) flavors, the /x^^'^'^^ are RGE- 
dependent coefficients, and Xi (X2) are the up- type (down-type) Yukawa matrices. 
The matrices Bi = U*U^, appearing in the equations in Sec. || above, are obtained 
from the f/j matrices that diagonalize Xfj^, and the f/j matrices that diagonalize the 
right-handed quark mass matrices. Because of the simple form for Xfj^ in Eq. (^Dp, 
it can be shown that f/j = f/j and therefore Bi = I 0. [Note that the quark mixing 
matrices f/j and Vi mentioned in this section are different from the chargino mixing 
matrices Uij and Vij mentioned in Section 2.] 

The other relevant set of matrices are F, = V^V^j , obtained from the Vi matrices 
that diagonalize Xfj^ and the Vi matrices that diagonalize the left-handed quark mass 
matrices. In the case of A^ ^ A;,, which requires tan/3 ~ 1, the X2X2 oc A^ term 

^The sign of these corrections depends on the observable being calculated: in the case of B{t -^ 
cV) they are positive, whereas in one-loop supersymmetric corrections to <j{pp — > tt) pO[ they were 
negative. 



in Eqs. (|51|j5^ ) is small compared to the XiXl oc A^ term, and therefore the former 
may be neglected. This implies that both X^^ and X21 are diagonalized by the same 
matrix: V^ = V2 = Vi, and therefore Fi = 14 ^1 = ^ whereas r2 = V2V2 = V1V2 = K 
reduces to the regular CKM matrix 0. As the quark-squark-neutralino couplings in 
flavor space are proportional to Fj, we see that for Aj ^ A;, there are no flavor off- 
diagonal couplings in the up-quark sector, as required for a unsuppressed contribution 
to the 'neutral' electroweak-like contributions to t ^ cV. 

One might consider instead a scenario where At ~ Xf,, as would be consistent 
with tan/3 ^1. In this case the X2X2 oc A^ term in Eqs. (|5T|j5^ is no longer 



negligible and Fi 7^ / is expected. The precise form of Fi requires a complicated 
calculation, essentially solving the matrix renormalization group equations. For our 
purposes here it suffices to consider the following effective form 



Fi = 1 e , (53) 




where e parametrizes the size of the ratio \b/\t- For moderate values of tan/5 this 
form should be adequate {i.e. e not too close to 1). We still expect F2 ~ K. We 
assume that the lower (2 x 2) right corner of Vi is approximately the identity to relate 
Tql to Fi. 

The above forms for Fi 2 plus the result B12 = I above, allow us to evaluate 
numerically the branching ratios of Sec. |^. Perhaps the most optimistic top factory 
being contemplated at the moment is a high-luminosity upgrade of the Tevatron, 
where studies show that one might be sensitive to B{t —>■ cy) ~ 4x 10~^ (8x 10~^) |14] 



B{t ^ cZ)^Ax 10^3 (6 X 10"^) 13, and B{t ^ eg) ^ 5 x IQ-^ (1 x IQ-^^) with 
an integrated luminosity of 10 (100) fb~^, where the branching ratios are with respect 
to F(t — > bW). These expected sensitivities will not allow direct tests of the Standard 
Model predictions for these processes: B{t -^ 07)^^ ~ 10"^^ B{t -^ cZf^ ~ 10"^^ 
and B{t — > cg)^^ ~ 10"^'' P, but might uncover virtual new physics effects that 
enhance these rates over Standard Model expectations. 

Indeed, we generally find that B{t ^ cV) greatly exceeds the corresponding 
Standard Model contribution, but unfortunately falls below the expected experimen- 
tal sensitivities, as was observed also in previous studies of the QCD-like corrections 
IQ. Specifically, concentrating on the 'neutral' electroweak-like corrections, we have 
as the dominant inputs the masses of the charm-squark and top-squark, the top mix- 
ing angle, the mass of the neutralino(s), and the neutralino composition. The results 
for V = g,'y scale with e^ as defined in Eq. (|53|) ; we take e = 0.5 for concreteness. (For 
V = Z, the cross section increases montonically with e for e < 0.5.) Numerically! we 
find that when rrit ~ m^ + mi^ , the branching ratios are enhanced compared to off- 
resonance values by a factor of 2-10. This factor depends on the specific combination 
of neutralino and top-squark masses that satisfy this relation (all other parameters 

^We used the software package f f llq] to evaluate the Passarino-Vehman functions. 



being kept fixed); the enhancement decreases with increasing stop mass and so is 
maximised when m^^ is 60 GeV. The off-resonance values themselves are larger than 
the Standard Model predictions for not-too-heavy sparticles. We also verify that 
large m^ — m^^ mass splitting enhances the results, because of its GIM-violating ef- 
fect. (A useful test of our code is that the branching ratios go to zero due to the GIM 
mechanism, if we set squark masses equal.) With regards to neutralino composition, 
the largest branching ratios are obtained for neutralinos with comparable bino and 
higgsino admixtures. Increasing the scale of the sparticle masses typically leads to a 
rapid decrease in the branching ratios for g and 7. For the Z the decrease is very 
gradual, as has also been noted in Ref. f^. For g and 7 the cross-section seems to 
be maximised for stop mixing angles close to or vr. The effects of mixing for Z are 
very dependent on the other parameters chosen, such as the neutralino composition, 
etc. The effect of varying tan/5 is of 0(1). 

In varying the different parameters above we have worked in the most general 
framework of the MSSM, in which the various parameters can be varied independently. 
In a more specific model, such as one with universal scalar masses and radiative 
electroweak symmetry breaking, these parameters are not all independent. Although 
our choice of mixing matrices was motivated by certain specific scenarios we vary 
our parameters freely so as to look for the maximal supersymmetric contributions. 
Furthermore, we choose mc^ 2 and rrif^ to be ~ ITeV. 

For the most optimistic values of the parameters, i.e., when the above enhanc- 
ing circumstances all simultaneously occur, we find B{t — > 07) < 2 x 10~'^, B{t — > 
cZ) <4 X 10"'^ and B{t ^ cg)<3 x 10"^. We see that B(t -^ eg) is the one closest 
to the level of experimental sensitivity expected at the Tevatron, so perhaps it would 
be the mode to be first observed at a future sufficiently sensitive machine. This hope 
is further enlarged by recalling that i?(t — > eg) receives comparable contributions 
from the QCD-like supersymmetric corrections @], which we have not evaluated here. 
Eventually, such a process can be a possible test for supersymmetry. 

We have also evaluated the charged electroweak-like corrections and have found 
them to be smaller (typically by a factor of 10 or more) than the neutral electroweak- 
like corrections discussed above for g and 7, but a factor of 10 higher for the Z. 
(Again we assume that the lower (2 x 2) right corner of V2 is approximately the 
identity to relate Tql to r2.) In the case of universal squark masses at the unification 
scale, GIM-violating bottom-squark strange-squark mass differences are generated by 
RGE evolution, resulting in shifts to the left-handed down-like squark mass matrices. 
The dominant term is from the second term of Eq. ( ^2]) which may be rewritten 
as —\e\K'^{mu)'^K, where m^ = {mu,mc,mt} and \e\ < 1 is an RGE-dependent 
constant. Inserting the values of the CKM matrix elements we find (approximately): 
r?7,|^ —> 771?^ — I c|(rri( /5)^ and rri? —>■ m3 — |c|m^. Choosing the maximal (|c| = 1) mass 

splittings, we find B{t -^ 07) < 10-^ B{t ^ cZ) < 2 x 10-^ and B{t -^ eg) < 10"^ 
for squark masses as low as experimentally allowed. (These results are not much 
altered even if one drops the assumption of universal squark masses at the unification 
scale.) The numerical results for the charged electroweak-like corrections cannot be 



compared with the corresponding ones in Ref. [Q] because, as we explained above, the 
formulas presented in Ref. are inconsistent with gauge invariance constraints. 

We finally try to connect up with the recent literature in Ref. 0, where in 
addition to the t ^ eg decay mode, people have considered hadronic processes like 
pp —>■ tc, which might be more easily detectable. These works ignore any substructure 
that the t-c-g vertex might have, and replace it all by an effective scale A, defined for 
instance by the relation r(t -^ eg) = Sasfn^/SA'^. A given branching ratio obtained 
in the supersymmetric theory then corresponds to a scale A in the effective theory. 



Dividing this expression by r(t -^ bW) we find that A ^ 1 TeV/JB{t -^ eg). There- 
fore a supersymmetric prediction of B(t -^ eg) ~ 10^^ corresponds to A ~ 300 TeV. 
The point of this exercise is to note how misleading such estimates of new-physics 
scales might be, as this actually corresponds in our case to sparticle masses of a few 
hundred GeV! 
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Figure 1: Feynman diagrams for one- loop electroweak-like supersymmetric contribu- 
tions to the t-c-V {y = 7, Z, g) vertex. In the figure x represents the chargino or 
the neutrahno and q represents the down-type or up-type squarks respectively. The 
subscripts are explained in the text. The arrows on the squark lines indicate the 
direction of fiow of flavor; the arrows on the gauge bosons indicate the direction of 
momentum flow. Diagram (b) is absent for V = g, and for V = '~f when x = X^- 
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